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SUMMARY 


Hydrogen  ingress  in  three  high-strength  alloys,  4340  steel.  Monel  K500, 
and  MP35N,  was  shown  to  occur  in  all  cases  under  interface  control;  under  the 
test  conditions  used,  the  rate  of  ingress  was  controlled  by  the  flux  across 
the  alloy  surface  into  the  bulk  and  not  by  hydrogen  diffusion  in  the  alloy. 
Under  interface  control  conditions,  the  hydrogen  ingress  model  was  valid  for  a 
metal  independent  of  the  electrolyte  used,  provided  that  the  electrode  surface 
is  unaffected  by  (1)  the  electrolyte  due  to  corrosion  or  a  change  in  any  film 
present,  or  (2)  cathodic  charging  due  to  bubble  formation  or  film  reduction. 

The  apparent  rate  constant  (ka)  for  hydrogen  trapping  in  4340  steel  was 
determined  for  a  range  of  electrolytes,  but  the  values  of  k  =  obtained  in  all 
cases  except  an  acetate  buffer  (pH  4.8)  were  affected  by  reactions  involving 
the  metal  surface  itself.  The  irreversible  trapping  constant  (k)  for  the 
steel  was  therefore  obtained  from  kfl  for  the  acetate  and  was  evaluated  by 
estimating  the  effect  of  reversible  traps  from  hydrogen  diffusivity  data.  The 
value  of  k  for  the  4340  steel  was  4  s-1,  and  the  density  of  irreversible  traps 
derived  from  k,  2  x  10®  m~®,  indicated  that  the  primary  irreversible  traps 
were  MnS  inclusions.  The  same  value  of  k  and  therefore  trap  density  was 
obtained  for  the  steel  irrespective  of  heat-treatment,  further  indicating  that 
the  irreversible  traps  are  stable  entities  such  as  MnS  inclusions. 

The  apparent  trapping  constant  for  Monel  K500  was  independent  of  the 
electrolyte  and  was  used  with  diffusivity  data  to  determine  the  irreversible 
trapping  constant.  The  value  of  k  was  0.042  s-1  which  gave  a  trap  density 
three  orders  of  magnitude  less  than  the  trap  density  expected  on  the  basis  of 
the  concentration  of  sulfur  atoms.  The  value  of  the  calculated  trap  density 
was  attributed  primarily  to  sulfur  segregated  as  clusters  at  grain  boundaries 
and  secondarily  to  the  effects  of  the  surface  film,  which  may  well  have  a 
different  sulfur  concentration  than  the  alloy. 
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Unlike  that  for  the  Monel,  the  trapping  constant  measured  for  MP35N 
appeared  to  increase  with  interfacial  pH.  This  increase  is  thought  to  be 
related  to  a  change  in  film  thickness  with  increasing  pH  of  electrolyte. 
Because  of  a  lack  of  data  for  the  diffusivity  of  hydrogen  in  MP35N,  the  rate 
constant  for  irreversible  trapping  was  approximated  to  the  measured  value.  As 
for  the  Monel,  the  trap  density  calculated  from  k  was  three  orders  of 
magnitude  less  than  expected  on  the  basis  of  the  sulfur  and  phosphorus  content 
in  the  alloy.  This  difference  was  similarly  ascribed  primarily  to  sulfur  and 
phosphorus  segregated  at  the  grain  boundaries  and  secondarily  to  effects  of 
the  surface  film. 


The  two  nickel-concaining  alloys,  Monel  K500  and  MP35N,  exhibit  a  lower 
flux  of  hydrogen  across  the  metal  surface  and  a  lower  irreversible  trapping 
constant  than  the  steel.  The  difference  in  both  the  nature  of  the  principal 
irreversible  traps,  which  were  determined  from  the  trapping  constants,  and  the 
interfacial  flux  can  be  expected  to  contribute  substantially  to  the  difference 
in  susceptibility  of  the  steel  and  the  two  nickel-based  alloys  to  hydrogen 
embrittlement . 
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INTRODUCTION 


The  ingress  of  hydrogen  into  steel  and  other  alloys  undergoing  corrosion 
can  severely  affect  the  properties  of  these  alloys,  especially  high-strength 
steels.  In  contrast,  certain  high-strength  nickel-based  alloys,  such  as  Monel 
K500  and  MP35N,  appear  to  be  much  more  resistant  to  hydrogen  embrittlement, 
although  some  failures  of  these  two  alloys  because  of  hydrogen  embrittlement 

1  9 

under  severe  conditions  have  been  reported.  ' 

Structural  features  such  as  dislocations,  internal  interfaces,  and  grain 
boundaries  are  potential  trapping  sites  for  diffusing  hydrogen.  The  hydrogen 
atoms  are  partitioned  between  the  traps  and  normal  lattice  sites  depending  on 
the  energy  binding  the  hydrogen  atom  to  the  trapping  site,  and  traps  are 
usually  classified  as  reversible  or  irreversible.  The  nature  of  the  traps 
within  a  metal  critically  influences  the  susceptibility  of  the  metal  to 
hydrogen  embrittlement. 

A  model  was  developed  previously3  to  represent  the  diffusion  and  trapping 
of  hydrogen  atoms.  The  model  is  coupled  to  the  use  of  a  potentiostatic 
double-pulse  (PDP)  technique  that  has  been  applied4  to  iron  in  H2S  solutions 
and  was  shown  to  provide  quantitative  information  on  both  the  mobile  and 
trapped  hydrogen  in  the  metal.  The  technique  is  a  transient  type  rather  than 
one  involving  steady-state  measurements  more  commonly  used  for  electrochemical 
studies  of  hydrogen  ingress.  Only  one  other  group3  has  reported  use  of  this 
technique,  and  their  data  analysis  did  not  include  hydrogen  trapping. 

In  this  study,  the  technique  and  the  associated  model  were  applied  to  a 
high-strength  steel  {AISI  4340)  and  two  high-strength  nickel-based  alloys 
(MP35N  and  Monel  K500)  to  determine  the  individual  hydrogen  ingress  and 
trapping  characteristics  and  thereby  obtain  some  insight  into  the  differing 
resistance  of  these  alloys  to  hydrogen  embrittlement. 
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Previous  electrochemical  investigations  of  hydrogen  ingress  have  been 
concerned  mainly  with  hydrogen  diffusion  through  thin  metal  membranes  using 


controlled  potential  or  current  conditions  at  both  interfaces.  In  contrast. 


the  PDP  technique  is  applied  to  sections  of  bulk  metal,  typically  but  not 


necessarily  in  rod  form.  The  metal  is  charged  with  hydrogen  at  a  constant 


potential  (E  )  for  a  time  (t  )  usually  in  the  range  0.5-40  s.  The  potential 


is  then  stepped  to  a  more  positive  value,  EA  (5-10  mV  negative  to  the 
open-circuit  potential) ,  resulting  in  an  anodic  current  transient  with  a 


charge  qa  due  to  reoxidation  of  H  atoms  at  the  same  surface,  according  to  the 


reaction 


metal 


— >  Hads - >  H  +  e 


where  HacJs  represents  hydrogen  atoms  adsorbed  at  the  metal  surface. 


The  effect  of  trapping  on  the  diffusion  of  hydrogen  is  taken  into  account 
by  modifying  Fick's  second  law  of  diffusion.  The  diffusion/trapping  model^ 
used  with  the  PDP  technique  was  based  on  a  single  trapping  term,  kc,  in  which 
the  rate  of  trapping  is  assumed  to  be  proportional  to  the  concentration. 


c(x,t),  of  diffusing  H  at  point  x  and  time  t: 


dc/dt  =  DL(d^c/dx^)  -  kc 


where  DL  is  the  lattice  diffusivity  of  hydrogen.  In  the  simplest  case,  it  is 
assumed  that  the  trapping  is  entirely  irreversible,  so  there  is  no  significant 


release  of  H  atoms  from  the  traps.  However,  equation  (2)  can  be  expanded  to 


allow  for  the  effect  of  reversible  traps  on  the  trapping  rate  constant  (k) . 
lino®  represented  the  combination  of  irreversible  and  reversible  trapping  in 
the  diffusion  equation  as  follows: 


>  i>kT»k* 


ms®, 


(3) 


dc/dt  +  N.(dn-/dt)  +  N  (dn  /dt)  =  D,(d2c/dx2) 
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where 


(dnr/dt)  =  krc(l-nr)  -  pnr 


(4) 


and 


(dn^/dt)  =  k^cd-n^)  (5) 

The  parameters  kr  and  p  are  the  trapping  and  release  rate  constants  for 
reversible  traps,  N  is  the  density  of  traps,  n  is  the  fraction  of  traps 
occupied,  and  the  subscripts  i  and  r  refer  to  irreversible  and  reversible 
trapping,  respectively.  In  his  treatment,  lino  assumed  a  low  trap  coverage  in 
both  the  irreversible  and  reversible  cases,  so  n^«l  and  nr«l.  Therefore, 
equation  (4)  reduces  to  (dnr/dt)  =  k^c  -  pnr  and  equation  (5)  to 
(dni/dt)  =  k^c . 

n 

It  is  further  assumed'  that  there  is  local  equilibrium  between  the 
trapped  hydrogen  and  the  lattice  hydrogen  at  the  reversible  sites;  that  is, 
(dnr/dt)  =  0  for  a  given  value  of  c.  Therefore,  nr  =  krc/p,  and  because 
c  =  f(t),  the  rate  of  reversible  trapping  is  given  by 

(dnr/dt)  =  kr/p(dc/dt)  (6) 

The  diffusion  equation  (3)  can  then  be  expressed  as 

dc/dt  +  N^k^c  +  (dc/dt) Nrkr/p  =  DL(d2c/dx2)  (7) 

and  therefore 


(dc/dt) (1  +  Nrkr/p)  +  Nikic  =  DL(d2c/dx2) 


(8) 
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McNabb  and  Foster"  derived  an  apparent  diffusion  coefficient  given  by 
Dl/ (1  +  Nrkr/p) ,  and  it  is  evident  that  division  of  equation  (8)  by 
(1  +  Nrkr/p)  will  similarly  give  an  apparent  diffusion  coefficient,  D  ,  and 
also  yield  an  apparent  or  measured  trapping  constant,  k, : 

a 

dc/dt  +  NjkjC/d  +  Nrkr/p)  =  [Dl/(1  +  Nrkr/p)  ]  (d2c/dx2)  ( 


dc/dt  +  kflc  =  Da(d2c/dx2) 


where  ka  =  k/ (1  +  Kr)  and  Da  =  D^/ (1  +  Kr) .  The  irreversible  trapping 

constant,  k,  is  given  by  N^k^,  and  Kr  is  an  equilibrium  constant  defined  as 

Nrkr/p.  Equation  (10)  has  also  been  derived  elsewhere  by  a  similar 
q 

treatment . 

The  effect  of  reversible  traps  on  k  can  therefore  be  taken  into  account 
through  the  equilibrium  constant  for  reversible  traps.  Because  the 
equilibrium  constant  also  relates  the  apparent  diffusivity  to  the  lattice 
diffusivity,  the  value  of  Kr  can  be  determined  from  diffusivity  data  if  such 
data  are  available  for  the  alloy  of  interest  and  for  the  lattice  of  the 
primary  components  of  the  alloy. 


On  the  basis  of  the  relatively  short  charging  times  (<40  s)  used  in  the 
PDP  technique,  the  irreversible  trap3  are  believed  not  to  become  saturated; 
that  is,  the  experimentally-determined  irreversible  trapping  constant,  k,  does 
not  vary  with  charging  time.  The  expression  for  the  rate  of  trapping 

O  » 

postulated  by  McNabb  and  Foster,"  k  N(1  -  n)c,  is  consistent  with  equation  (2) 
if  it  is  assumed  that  (1  -  n)  is  constant  on  the  time  scale  of  the 
potential-step  experiments.  The  irreversible  trapping  constant  is  then  given 

I 

by  k  =  k  N(1  -  n) ,  which,  for  low  trap  coverage,  corresponds  to  the  expression 
Ni^i  <?iven  above. 

Furthermore,  the  reversible  traps  are  assumed  to  release  all  their 
hydrogen  during  the  anodic  transient,  t  ,  so  that  all  the  hydrogen  remaining 


in  the  metal  at  the  end  of  ta  is  truly  irreversibly  trapped.  In  the  case  of 
4340  steel,  the  assumption  of  highly  reversible  traps  is  based  on  a  value  of 
500  estimated  for  Kr  (see  below) .  The  density  of  reversible  traps  in  4340 
steel10  is  typically  in  the  range  lO^-lO27  m-2 .  Even  at  the  lower  limit  of 
this  range,  kr/p  =  5  x  10  m  ,  which  implies  that  the  release  constant  is 
large,  and  so  the  release  rate  will  be  high  even  when  the  fraction  of  traps 
occupied  is  very  low.  In  the  case  of  the  two  nickel-based  alloys,  the  binding 
energy  for  hydrogen  to  defects  such  as  vacancies  or  edge  dislocations  in  an 
fee  lattice  is  about  0.1  eV;  in  comparison,  the  activation  energy  for 
diffusion  in  this  type  of  lattice  is  roughly  0.4  eV.11  Therefore,  transport 
of  hydrogen  in  these  alloys  is  limited  by  diffusion,  with  little  hindrance 
from  reversible  traps. 

The  form  of  equation  (10)  is  identical  to  that  of  equation  (2),  which  has 
been  solved  analytically2  for  two  cases  characterized  by  the  rate-determining 
step  for  the  transfer  of  Hads  into  the  metal  during  the  charging  time.  In 
both  cases,  it  is  assumed  that  that  the  layer  of  Hads  adjusts  very  rapidly  to 
a  potential  step.  For  a  given  value  of  the  charging  potential,  Ec,  the 
hydrogen  ingress  flux  into  the  metal  is  considered  constant.  If  equilibrium 
is  established  between  the  concentration  of  hydrogen  in  the  metal  at  x  =  0  and 
the  adsorbed  hydrogen  layer,  the  ingress  of  hydrogen  will  be  controlled  by 
diffusion  in  the  bulk  metal.  However,  if  the  flux  across  the  interface  is 
very  slow,  the  rate  of  ingress  will  be  determined  by  the  flux.  In  this  case, 
hydrogen  ingress  occurs  under  interface  control. 
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EXPERIMENTAL  PROCEDURE 


The  chemical  composition  of  each  alloy  was  provided  by  the  manufacturer 
and  is  given  in  Table  1 . 

Table  1 

CHEMICAL  COMPOSITION  OF  THE  THREE  ALLOYS  TESTED 


Element  4340  steel  Monel  K500 


MP35N 


<0.01 


0.009 


20.19 


0.031 


0.001 


Bal  =  balance,  n  =  not  present  or  not  given. 


The  4340  steel  was  supplied  as  1.27-cm  plate  produced  from  an 
argon-oxygen-decarburized  heat  and  subsequently  vacuum-arc  remelted  (VAR) 
before  being  forged,  rolled,  and  annealed.  The  as-received  plate  was  machined 
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into  two  rods  1.25  cm  in  diameter  and  5  cm  long,  which  were  heat  treated  to 
Rockwell  C  hardness  (HRC)  values  of  41  and  53.  These  values  correspond  to 
yield  strengths  of  175  and  260  ksi,  respectively.  The  conditions  of  heat 
treatment  are  shown  in  Table  2. 


Table  2 

CONDITIONS  OF  HEAT  TREATMENT  FOR  4340  STEEL 


Step 

HRC  41 

HRC  53 

Austenized 

830°C  for  45  min 

850°C  for  45 

min 

Quench 

Oil 

Oil 

Tempered 

450°C  for  60  min 

180°C  for  80 

min 

Quench 

Water 

Water 

10-1  ? 

MnS  inclusions  are  well-recognized  as  strong  hydrogen  traps,  and  a 

thermal  analysis  study15  of  4340  steel  showed  that  the  temperature  for 
hydrogen  release  increases  in  the  following  order:  ferrite-carbide  interface  < 
dislocation  <  microvoid  <  MnS  interface.  Accordingly,  MnS  inclusions  are 
likely  to  be  principal  irreversible  trapping  sites  in  4340  steel.  In  an 
earlier  study  in  which  the  nature  and  distribution  of  inclusions  in  the  4340 
steel  used  in  this  work  were  determined, 15  it  was  found  that  there  were 
~105  itT^  MnS  inclusions  randomly  distributed  through  the  material.  The 
inclusions  were  ellipsoidal  with  the  major  axis  transverse  to  the  rolling 
direction.  On  the  basis  that  the  inclusions  were  spherical,  the  concentration 

Q 

of  inclusions  was  calculated  to  be  2  x  10  m  and  the  mean  radius  of  the 
inclusions  averaged  over  the  transverse  and  parallel  directions  was  9.2  (im. 

The  Monel  K500  was  supplied  by  Huntington  Alloys  as  an  unaged,  1.27-cm, 
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cold-drawn  rod.  The  HRC  value  of  the  as-received  (AR)  material  was  26,  which 
corresponds  to  a  yield  strength  of  -110  ksi.  A  section  of  the  rod  was 
age-hardened  to  HRC  35,  giving  a  yield  strength  of  about  159  ksi.  The 
age-hardening  was  performed  according  to  the  procedure  recommended  by  the 
manufacturer:  Heat  for  8  h  at  600°C  and  then  cool  to  480°C  at  a  rate  of  ll°/h. 


MP35N  is  based  on  the  quaternary  of  cobalt,  nickel,  chromium,  and 
molybdenum.  It  is  a  vacuum-induction,  vacuum-arc  remelted  alloy  and  was 
supplied  by  Latrobe  Steel  Co.  as  a  1.36-cm  cold-drawn  rod  that  had  been  aged 
by  heating  at  590°C  for  4  h  and  then  air-cooling.  The  hardness  and  yield 
strength  for  the  as-received  material  were  HRC  45  and  269  ksi,  respectively. 


The  test  electrodes  of  each  alloy  were  fabricated  from  5-cm  lengths  of 
rod.  A  stainless  steel  rod  3  mm  in  diameter  screwed  into  one  end  of  the  rods 
to  support  the  electrode  in  the  cell  and  provide  electrical  contact.  Each  rod 
was  press-fitted  into  a  Teflon  sheath  so  that  only  the  planar  end  surface  of 
the  rod  was  exposed  to  the  electrolyte.  The  surface  was  polished  with 
successively  finer  grades  of  SiC  paper  from  320  to  600  grit  followed  by 
0.05-(lm  alumina  powder. 


Cell 

A  conventional  three-electrode  electrochemical  cell  was  used.  The 
counter  electrode  was  a  graphite  rod,  and  the  reference  electrode  assembly 
consisted  of  a  saturated  calomel  electrode  (SCE)  inserted  in  a  Luggin  probe. 
All  potentials  are  given  with  respect  to  the  SCE.  Various  electrolytes 
ranging  in  pH  from  2.9  to  9.0  (Table  3)  were  used. 


Table  3 

COMPOSITION  AND  pH  OF  TEST  ELECTROLYTES 


Designation 


Composition 


Sulfuric  acid 

0.0009  mol  L_i  H2S( 

Acidified  chloride 

3  wt%  NaCl 

3.7  x  10"5  mol  L-1 

Acetate  buffer 

1  mol  L-1  CH3COOH 

1  mol  L"1  CH3COONa 

Phosphate  buffer 

0.5  mol  L-1  Na2HPO 
0.5  mol  L-1  NaH2PO 

Bicarbonate  buffer 

0.5  mol  L-1  NaHC03 
0.01  mol  L-1  NaOH 

Measured  values 


Each  electrolyte  contained  15  ppm  As2C>3  as  a  hydrogen  entry  promoter.  The 
electrolyte  was  deaerated  with  argon  for  1  hour  before  measurements  began  and 
throughout  data  acquisition.  All  measurements  were  performed  at  22  ±  2°C. 


The  PDP  technique  was  performed  using  a  Princeton  Applied  Research  Model 
173  potentiostat  with  a  Model  276  interface  coupled  to  a  Macintosh 
microcomputer  via  an  Iotech  Mac  488A  bus  controller.  A  BASIC  program, 

PDP  173/488  (Appendix  A),  was  developed  to  provide  the  perturbation 
(double-pulse)  waveform,  set  control  parameters,  and  acquire  current  data.  Th 
open-circuit  potential  (EQC)  of  the  test  electrode  was  sampled  before  each 
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charging  time  (tc)  and  the  electrode  was  then  potentiostatically  controlled  at 
Ea,  which  was  set  to  10  mV  negative  of  EQC.  This  approach  allows  for  any 
drift  in  EQC  and  therefore  reflects  the  corrosion  state  of  the  metal  at  all 
times  during  the  range  of  EQ  and  tc  values  used.  The  practical  consequence  of 
this  procedure  is  that  Eft  is  always  maintained  slightly  negative  of  the 
existing  EQC,  and  thus  anodic  oxidation  of  the  metal  is  minimized. 

Control  Parameters 

The  principal  control  parameters  are  the  charging  potential  (Ec) ,  the 
period  of  the  anodic  and  cathodic  transients  (ta  and  tc,  respectively) ,  and 
the  current  sampling  rate.  These  parameters  have  several  electrochemical  and 
metallurgical  limitations. 

The  primary  concern  is  to  generate  a  sufficiently  rapid  discharge  of  H+ 
during  t  to  transfer  a  significant  quantity  of  hydrogen  into  the  metal. 
However,  it  is  undesirable  for  molecular  hydrogen  to  form  to  an  extent  that  an 
appreciable  number  of  bubbles  develop  on  the  surface,  because  the  bubbles 
would  influence  the  surface  coverage  of  Hadg  and  also  reduce  access  of  the 
electrolyte  to  the  surface.  These  factors  are  determined  by  the  hydrogen 
evolution  kinetics  of  the  alloy,  and  consequently  the  magnitude  of  the 
cathodic  potential  step  was  varied  according  to  the  alloy  being  tested. 

The  cathodic  charging  times  were  selected  to  vary  from  0.5  to  40  s,  which 
provided  a  range  covering  almost  two  decades.  With  charging  times  under 
0.5  s,  insufficient  hydrogen  enters  the  metal,  and  with  very  short  charging 
times,  the  current  is  dominated  by  charging  of  the  double  layer  at  the 
surface.  With  prolonged  charging  of  the  metal,  irreversible  traps  may  become 
saturated,  and  the  mathematical  model  used  in  the  analysis  would  be 
invalidated.  Moreover,  excessively  long  charging  times  might  cause  internal 
damage  to  the  metal  due  to  the  accumulation  of  hydrogen.  The  reproducibility 
of  anodic  charge  data  in  consecutive  tests  with  tc  =  40  s  suggested  that  traps 
were  neither  saturated  nor  generated  as  a  result  of  structural  damage. 
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The  time  of  data  acquisition  (ta)  during  the  anodic  transient  was  made 
equal  to  that  for  the  cathodic  transient  on  the  basis  that  this  time  would  be 
long  enough  for  the  deepest  untrapped  hydrogen  to  diffuse  back  to  the  metal 
surface  and  for  hydrogen  in  highly  reversible  traps  to  be  released,  as 
discussed  above.  Accordingly,  it  was  assumed  that  qa  is  equal  to  the  anodic 
charge  q'  (°°)  corresponding  to  ta  =  00  (see  below)  . 

The  sampling  time  for  the  current  was  usually  selected  to  allow 
sufficient  resolution  of  the  current  transient.  Various  sampling  times  were 
tested.  Table  4  shows  the  values  selected  and  the  total  number  of  data  points 
acquired  during  the  anodic  and  cathodic  transients  for  each  charging  time;  the 
anodic  and  cathodic  parts  each  contain  half  the  number  of  points  shown. 


Table  4 

DATA  ACQUISITION  CHARACTERISTICS 


Sampling  time  (ms)  Number  of  points 


RESULTS 


The  anodic  current  transients  for  all  three  alloys  characteristically 
exhibited  the  same  shape  over  the  range  of  charging  potentials  and  times  of 
interest.  A  typical  current/time  (i/t)  transient  is  shown  in  Fig.  1.  In  all 
cases,  analysis  of  the  data  in  terms  of  the  diffusion-control  model  indicated 
that  (1)  ka  varied  according  to  the  number  of  points  analyzed  from  the 
transient  for  each  value  of  tc,  (2)  ka  decreased  with  increasing  tc  for  a 
given  potential,  and  (3)  kg  generally  decreased  with  an  increase  of  potential 
for  a  given  tc.  The  diffusion  control  model  therefore  was  not  valid  for  any 
of  the  three  alloys,  but  in  all  cases  the  data  could  be  analyzed  in  terms  of 

the  interface  control  model.  According  to  the  interface  control  model,  the 
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total  charge  passed  out  is  given  m  nondimensional  form  by 


q'(co)  =  Vr{1  -  e-R/V(7tR)  -  [1  -  1/ (2R)  ]erf^R}  (11) 

The  nondimensional  terms  are  defined  by  Q*  (°°)  =  q’ (°°)  /  [FjV  (tc/ka)  ]  and 

R  =  katc  where  J  is  the  ingress  flux  in  mol  cm  s  and  the  anodic  charge 

.  —7 

q'  (°°)  is  equated  to  qa  in  mC  cm  .  For  the  model  to  be  applicable,  it  must  be 
possible  to  determine  a  trapping  constant  for  which  the  ingress  flux  is 
constant  over  the  range  of  charging  times.  The  data  were  analyzed  using  the 
BASIC  program,  INTERFACE,  given  in  Appendix  B.  The  values  for  the  trapping 
constant  and  flux  can  be  used  to  calculate  the  amount  of  charge  (qT) 
irreversibly  trapped  under  conditions  of  interface  control  for  which  the 
trapped  charge  is  given  nondimensionally  by 


Qt  =  [Vr  -  1/ (2^R) ]erf^R  +  e-R/Vm 


The  data  analysis  using  the  interface  control  model  was  restricted  to 
charging  times  of  10  to  40  s.  There  are  two  reasons  for  this  restriction: 
(1)  At  short  charging  times,  the  values  of  qa  may  include  a  significant 
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contribution  due  to  oxidation  of  the  adsorbed  hydrogen  layer,  qads,  and  the 
model  assumes  that  qa  is  associated  entirely  with  oxidation  of  hydrogen  that 
entered  the  metal.  (2)  The  cathodic  current  during  tc  changes  to  some  extent 
in  its  early  stages,  and  this  change  may  affect  the  flux  of  hydrogen  across 
the  metal  surface.  The  model  assumes  that  mass  transport  of  H+  to  the  metal 
surface  during  t  is  constant,  and  this  assumption  is  more  justified  for 
values  of  tc  510  s,  where  the  current  is  approximately  constant  for  mo s t  of 
the  charging  time. 


The  anodic  charge  as  a  function  of  t  could  generally  be  analyzed  with 
the  interface  control  model  to  the  extent  that  values  of  the  apparent  trapping 
constant  and  flux  were  obtainable.  However,  the  steel  was  susceptible  to  film 
formation  or  metal  dissolution  in  all  the  electrolytes  except  the  acetate 
buffer.  In  such  cases,  the  validity  of  the  results  was  uncertain  because 
anodic  reactions  involving  the  alloy  itself  appeared  to  introduce  significant 
errors  into  the  values  of  qa  for  H  atom  oxidation,  and  thereby  cause  anomalies 
in  the  results  for  Jca  and  J. 

Consistent  values  were  obtained  with  the  acetate  buffer,  which  proved  to 
be  the  most  suitable  electrolyte  for  this  type  of  study  because  it  provided  an 
adequate  and  steady  supply  of  H+  for  reduction  but  was  not  so  acidic  as  to 
cause  substantial  attack  of  the  alloy.  Accordingly,  the  results  for  this 
electrolyte  are  presented  first  for  all  three  alloys  tested. 

Acetate .  Values  of  qa  were  plotted  against  Vtc,  which  was  used  to  be 
consistent  with  the  previous  work,^'^  and  data  for  a  potential  step, 

T\  (Ec  -  Eoc) ,  of  -0.4  V  [Ec  =  -0.985  V  (SCE) ]  are  shown  in  Fig.  2.  The  anodic 
charge  (qads)  arising  from  oxidation  of  the  adsorbed  H  layer  can  be  estimated 
by  extrapolation  to  tQ  =  0,  and  therefore,  data  are  displayed  only  for  shorter 
charging  times  (t  <15  s)  to  show  the  extrapolated  curve  more  clearly;  the  full 
range  of  data  is  presented  in  Fig.  3. 
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The  dependence  of  qa  on  tc  can  be  represented  by  a  second-order  polynomial: 


q,  =  7.9  x  10“5  +  0.031  t„1/2  +  0.011  t„  (13) 

|  ci  ^  C 

whereas  in  the  absence  of  trapping,  qg  is  expected  to  vary  linearly  with  tc. 
From  equation  (13),  it  is  apparent  that  qads  <10-4  mC  cm-2  which  can  be 
|  regarded  as  negligible.  Accordingly,  qa  is  assumed  to  be  associated  entirely 

i  with  absorbed  hydrogen. 

t 

i 

|  The  data  for  both  heat-treated  specimens  were  analyzed  using  the 

|  interface  control  model.  Values  of  the  interfacial  flux  and  apparent  trapping 

,  constant  determined  for  charging  times  from  10  to  40  s  are  given  for  each 

«  potential  step  in  Table  5.  In  both  cases,  current  transients  were  obtained  at 

i 

I  a  lower  overpotential  (-0.25  V),  but  the  results  for  qa  were  too  low  to 

1  determine  reliable  values  for  k  and  J. 


Table  5 


INTERFACIAL 

FLUX  AND 

TRAPPING  CONSTANTS  FOR 

IN  ACETATE 

4340  STEEL 

Type 

H  (V) 

-2  -1 

J  (mnol  cm  s  ) 

ka  (s'1) 

HRC  41 

-0.30 

0.045 

0.009 

-0.35 

0.07 

0.006 

0.19 

0.008 

-0.45 

0.36 

0.008 

HRC  53 

-0.30 

C  .  07 

0.008 

-0.35 

0.13 

0.007 

-0.40 

0.26 

0.008 

-0.45 

0.55 

0.009 

L*0 


For  both  heat-treated  specimens,  ka  is  essentially  independent  of  T|,  as 
expected,  because  the  density  of  traps  within  the  metal  should  not  be  affected 
by  the  electrochemical  surface  conditions.  Moreover,  in  both  cases  ka  has 
approximately  the  value  0.008  s-^.  Using  this  value  for  ka  and  the 
appropriate  value  of  J  in  equation  (11),  qa  was  calculated  for  charging  times 
from  0.5  to  40  s  and  compared  with  the  corresponding  experimental  results. 

The  two  sets  of  data  for  the  HRC  53  steel  charged  at  f|  =  -0.4  V  are  shown  in 
Fig.  3,  where  Q'  is  given  by  Q ' (°°) V (tc/k)  =  qa/FJ.  The  form  of  the  curves  and 
the  close  fit  between  them  are  typical  for  both  specimens.  The  agreement 
between  the  two  sets  of  data  indicates  the  independence  of  both  k  and  J  with 
respect  to  tc. 

In  contrast,  the  flux  does  increase  with  T|  as  a  result  of  the  dependence 
of  J  on  the  surface  coverage  of  Hacjs  and  therefore  on  potential.  The  surface 
coverage  and  hence  the  flux  should  exhibit  an  exponential  dependence  on 
potential,  which  is  in  fact  observed  for  both  specimens,  as  shown  in 
logarithmic  form  in  Fig.  4.  The  regression  coefficients  (r)  are  0.99  and 
1.00,  and  interestingly,  although  the  amount  of  data  is  limited,  the  two 
slopes  (m  =  -6.29  and  -5.97)  appear  to  be  close.  This  closeness  suggests  that 
the  rate  constant  for  the  flux  across  the  interface  is  determined  primarily  by 
electrochemical  characteristics  and  not  by  alloy  microstructure. 

The  dimensionalized  trapped  charge  (qT)  for  4340  steel  was  calculated 
using  equation  (12) .  The  data  for  the  HRC  53  specimen  charged  at  t\  =  -0.4  V 
are  given  in  Table  6,  which  also  shows  the  experimental  data  for  the  anodic 
charge  (qa)  and  cathodic  charge  (qc)  passed  during  the  respective  potential 
pulses  and  the  charge  associated  with  the  entry  of  hydrogen  into  the  metal 
(qin )  given  by  the  sum  of  qT  and  qa.  The  cathodic  charge  was  shown  to  be 
linearly  dependent  on  tc  (Fig.  5) .  This  linear  dependence  indicates  that  the 
cathodic  current,  and  therefore  supply  of  H+,  was  essentially  constant  even 
for  relatively  low  charging  times.  The  data  for  q^n,  qT,  and  qc  were  used  to 
obtain  two  ratios:  (1)  g-p^in'  corresponding  to  the  fraction  of  hydrogen  in 
the  metal  that  was  trapped;  and  (2)  gj_n/gc>  representing  the  fraction  of 
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charge  associated  with  hydrogen  entry  during  the  charging  step. 


Table  6 

CHARGE  PARAMETERS3  FOR  4340  STEEL  HRC  53  IN  ACETATE*3 


ft' 

1 

* 

y 

i 


i  0.5 

0.001 

0.028 

0.029 

0.021 

-0.62 

c 

0.046  S 

V 

1 

0.002 

0.044 

0.046 

0.037 

-1.23 

0.037  ft 

i 

,  2 

0.005 

0.064 

0.069 

0.069 

-2.37 

0.029  ft 

>  5 

i 

0.019 

0.123 

0.142 

0.132 

-5.67 

0.025  0 

10 

0.053 

0.212 

0.265 

0.199 

-10.96 

0.024  $ 

I  15 

0.096 

0.283 

0.379 

0.253 

-16.67 

0.023  | 

;  20 

0.146 

0.345 

0.491 

0.298 

-22.74 

0.022  5) 

;  25 

0.203 

0.441 

0.644 

0.315 

-29.03 

0.022  $ 

!  30 

0.265 

0.560 

0.825 

0.321 

-35.28 

0.023  5 

40 

0.401 

0.611 

1.012 

0.397 

-47.55 

0.021  | 

* 
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■  Monel  AR 

■  MP35N 


Despite  the  complex  variation  of  qT  with  tc,  has  an  approximately 

linear  dependence  on  Vtc,  as  shown  in  Fig.  7.  In  contrast,  the  ratio  <3j_n/qc 
is  independent  of  potential  and  tc  except  at  very  short  charging  times.  The 
fact  that  this  ratio  lies  between  0.02  and  0.03  (tc>l  s)  indicates  that  only 
about  2%-3%  of  the  hydrogen  atoms  formed  during  charging  enter  the  metal. 

H2.S&4-  The  electrode  surface  turned  grey  during  the  measurements, 
probably  as  a  result  of  etching  of  the  metal,  and  became  extensively  covered 
by  hydrogen  bubbles  even  at  relatively  low  overpotentials  (-0.15  V).  These 
two  factors,  especially  surface  blocking  by  the  bubbles,  were  assumed  to  be 
the  cause  of  meaningless  values  observed  for  qa.  The  anomalies  in  qa 
precluded  attempts  to  obtain  ka  and  J  for  the  steel  in  this  electrolyte. 

NaCl /HC1 .  In  this  electrolyte,  the  surfaces  of  both  steel  specimens 
exhibited  a  blue  tarnishing  characteristic  of  interference  effects  associated 
with  thin  films.  The  formation  of  a  film  may  have  been  due  to  the  increase  in 
pH  at  the  electrode  surface  because  of  the  unbuffered  nature  of  the 
electrolyte;  if  the  concentration  of  H+  at  the  electrode  surface  during 
reduction  became  low  enough,  the  pH  would  approach  neutral  and,  with  the 
reduction  of  H20,  even  become  alkaline. 

The  errors  in  qa  caused  by  the  formation  of  a  film  are  thought  to  be 
responsible  for  some  unexpected  features  in  the  behavior  of  ka .  Values  of  ka 
and  J  for  the  steel  in  acidified  NaCl  are  shown  in  Table  7,  and  it  is  evident 
that  ka  changes  erratically  with  potential,  especially  for  the  HRC  53 
specimen.  Moreover,  ka  for  the  HRC  53  steel  at  the  higher  overpotentials  is 
approximately  twice  as  large  as  that  for  the  HRC  41  specimen,  and  the  average 
value  of  ka  for  HRC  53  steel  is  0.042  s-^,  in  comparison  with  0.026  s-1  for 
the  lower  strength  steel.  Finally,  ka  for  both  specimens  is  high  in 
comparison  with  the  values  obtained  for  steel  in  the  acetate  buffer  of  the 
same  bulk  pH  as  the  NaCl/HCl  solution. 


Table  7 


INTERFACIAL 

FLUX  AND 

TRAPPING  CONSTANTS  FOR 

IN  NaCl/HCl 

4340  STEEL 

Type 

r\  (V) 

J 

_  I 

(nmol  cm  ^  s  ) 

ka  (a-1) 

HRC  41 

-0.25 

0.89 

0.036 

-0.30 

1.14 

0.028 

-0.35 

1.13 

0.019 

-0.40 

1.43 

0.029 

-0.50 

1.55 

0.019 

HRC  53 

-0.25 

0.44 

0.017 

-0.30 

0.75 

0.030 

-0.35 

1.48 

0.048 

-0.40 

2.08 

0.061 

-0.50 

3.92 

0.055 

Phosphate.  The  surface  of  the  steel  developed  a  dark  area,  which 
suggested  either  etching  or  the  formation  of  a  film,  and  although  trapping 
constants  could  be  obtained,  they  varied  randomly  with  increasing 
overpotential,  as  shown  in  Table  8.  As  with  the  NaCl/HCl  solution,  the 
average  values  of  ka  (0.043  s-'''  for  HRC  41  and  0.032  s-'1-  for  HRC  53)  are  high 
compared  with  those  for  the  acetate  buffer. 


Table  8 


INTERFACIAL 

FLUX  AND 

TRAPPING  CONSTANTS 

IN  PHOSPHATE 

FOR  4340  STEEL 

Type 

11  (V) 

J 

_  O  __  1 

(nmol  cm  s  ) 

ka  (s’1) 

HRC  41 

-0.25 

1.41 

0.057 

-0.30 

1.82 

0.020 

-0.35 

3.62 

0.046 

-0.40 

1.43 

0.049 

HRC  53 

-0.25 

0.70 

0.019 

-0.30 

1.90 

0.017 

-0.35 

5.74 

0.052 

-0.40 

7.51 

0.040 

Bicarbonate .  The  electrode  surface  generally  developed  species  and  often 
patches  of  blue  tarnish  typical  of  interference  films.  In  some  cases, 
tarnishing  was  not  visible  but  the  rate  constants  for  the  apparently 
untarnished  electrodes  fell  within  the  range  of  those  for  tarnished 
electrodes.  Moreover,  the  extent  of  visible  tarnishing  had  little  discernible 
effect  on  the  trapping  constant;  that  is,  the  results  were  within  the  limits 
of  experimental  error.  Table  9  shows  the  values  of  k=  and  J  obtained  for 

a 

electrodes  that  showed  a  typical  amount  of  tarnishing.  Again,  the  average 
values  of  ka  (0.028  s-1  for  HRC  41  and  0.050  s-1  for  HRC  53)  were  higher  than 
the  corresponding  results  obtained  with  acetate. 


Table  9 

INTERFACIAL  FLUX  AND  TRAPPING  CONSTANTS  FOR  4340  STEEL 

IN  BICARBONATE 


J  ( nmo 1  cm 


HRC  41 


HRC  53 


.50 

0.54  1 

.55 

0.79  1 

.60 

0.71  1 

.65 

0.86 

.50 

1.54 

.55 

1.52 

.35 

1.38 

i.40 

1.36 

0.048 

0.070 

0.045 


At  sufficiently  negative  overpotentials,  the  open-circuit  potential  of 
the  Monel  in  all  electrolytes  except  H2SO4  shifted  to  more  negative  values 
during  the  range  of  charging  times  used  and  extensive  bubble  formation 
occurred.  The  shift  in  EQC  probably  reflected  reduction  of  an  oxide  film  on 
the  alloy  surface;  the  implication  is  that  at  lower  values  of  T|,  the  reduction 
of  H+  and  the  subsequent  entry  of  H  atoms  into  the  metal  occurred  on  an 
oxide-covered  surface. 

Acetate .  The  anodic  charge  associated  with  oxidation  of  the  adsorbed 
layer  of  hydrogen  was  again  negligible;  in  this  case,  qac^s  =  7  x  10-<*  mC  cm~^ . 
Values  of  the  interfacial  flux  and  trapping  constant  determined  from  qa  data 
for  charging  times  from  10  to  40  s  are  given  in  Table  10. 
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INTERFACIAL  FLUX  AND  TRAPPING  CONSTANTS  FOR  MONEL  K500 

IN  ACETATE 


Type 

T)  (V) 

__o  _  1 

J  (nmol  cm  s  ) 

ka  (s"1) 

AR 

-0.65 

0.15 

0.027 

-0.70 

0.18 

0.027 

-0.75 

0.18 

0.019 

1 

o 

CD 

O 

0.21 

0.027 

HRC  35 

1 

O 

<T\ 

cn 

0.12 

0.016 

-0.70 

0.17 

0.018 

-0.75 

0.17 

0.012 

For  both  Monel  specimens,  the  results  exhibit  similar  patterns  to  those 
for  4340  steel:  (1)  ka  is  essentially  independent  of  Ec.  (2)  The  flux 
increases  with  T|,  but  not  to  the  same  extent  as  it  does  for  4340  steel.  Even 
at  the  high  overpotentials  used  for  the  Monel,  the  values  of  J  are  smaller 
than  those  for  the  steel,  presumably  because  of  the  effect  of  the  film. 


Values  of  Q'  calculated  using  the  appropriate  values  of  k.  and  J  for 

a 

charging  times  from  0.5  to  40  s  show  good  agreement  with  the  corresponding 
experimental  results,  as  shown  for  the  Monel  AR  and  T|  =  -0.8  V  [Ec  =  -0.905  V 
(SCE) ]  in  Fig.  8.  The  difference  between  the  two  sets  of  data  at  short 
charging  times  probably  reflects  the  influence  of  changes  in  the  cathodic 
current  in  its  early  stages  (see  below)  rather  than  a  dependence  of  k  and  J  on 
tc.  In  addition,  the  precision  of  the  qg  data  for  each  potential  was 
undoubtedly  affected  by  small  variations  of  up  to  15  mV  observed  in  EQC  over 
the  range  of  tc  and  particularly  by  the  more  pronounced  changes  in  EQC  at 


higher  overpotentials,  as  noted  above. 
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hydrogen . 

N aCl  /HC.l .  The  trapping  constants  for  the  HRC  35  specimen  showed  some 
scatter  over  the  range  of  overpotentials,  as  shown  in  Table  12,  but  the  mean 
value  of  0.029  s-'*'  was  consistent  with  that  for  the  acetate  solution. 


Table  12 

INTERFACIAL  FLUX  AND  TRAPPING  CONSTANTS  FOR  MONEL  K500  HRC  35 

IN  NaCl/HCl 


T|  (V) 

—2  —1 

J  (nmol  cm  s  ) 

ka  (s"1) 

-0.60 

0.16 

0.028 

-0.65 

0.23 

0.034 

1 

o 

o 

0.29 

0.040 

-0.75 

0.26 

0.019 

1 

o 

00 

o 

0.41 

0.016 

-0.85 

0.69 

0.023 

-0-90 

1.31 

0.057 

-0.95 

0.88 

0.015 

I 


Phosphate .  The  values  of  for  the  HRC  35  specimen  are  given  in 
Table  13.  They  show  somewhat  less  scatter  than  those  for  the  acidified  NaCl 
but  they  still  vary  from  0.010  to  0.023  s-^.  The  mean  value  is  0.017  s“  . 
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Table  13 

XNTERFACIAL  FLUX  AND  TRAPPING  CONSTANTS  FOR  MONEL  K500  HR C  35 

IN  PHOSPHATE 


J  (nmol  cm 


(s-1) 


.60 

0.19 

.65 

0.18 

.70 

0.22 

.75 

0.23 

.80 

0.23 

.85 

0.20 

.90 

0.27 

.95 

0.29 

0.014 
0.014 
0.018 
0.019 


The  values  of  ka  for  Monel  HRC  35  in  the  bicarbonate  buffer 


(Ta*\le  14)  are  similar  to  those  in  the  phosphate  buffer.  The  mean  value, 
_  -1 

0.019  s  ,  is  also  similar. 


FLUX 

Table  14 

AND  TRAPPING  CONSTANTS 

IN  BICARBONATE 

FOR  MONEL 

T\  (V) 

J  (nmol  cm-2  s-1) 

ka  (s"1) 

-0.60 

0.43 

0.024 

-0.65 

0.51 

0.023 

1 

o 

-J 

o 

0.50 

0.009 

-0.75 

0.89 

0.022 

-0.80 

0.85 

0.017 
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MP35N 

The  behavior  of  this  alloy  was  similar  to  but  more  marked  than  that  of 
the  Monel.  The  open-circuit  potential  for  MP35N  in  all  electrolytes  drifted 
considerably  over  the  range  of  tc  values  at  each  potential.  During  the  early 
stages  of  immersion,  EQC  showed  a  net  positive  drift  over  the  range  of 
charging  times  and  charging  potentials,  presumably  as  a  result  of  passivation 
of  the  freshly-polished  surface  of  the  electrode.  The  presence  of  a  film  was 
further  indicated  by  a  substantial  negative  shift  in  EQC  over  the  range  of 
charging  times  at  sufficiently  high  cathodic  overpotentials  and  by  an  increase 
in  the  values  of  qa  far  greater  for  these  overpotentials  than  for  lower 
overpotentials.  Thus,  as  with  the  Monel,  it  is  evident  that  hydrogen  enters 
MP35N  across  a  surface  film,  even  up  to  relatively  high  overpotentials. 

A  change  in  EQC  means  that,  for  a  given  T|,  Ec  is  not  constant  over  the 
range  of  charging  times.  Rather  than  sample  EQC  for  each  charging  time,  it 
could  have  been  sampled  only  at  the  start  of  each  new  overpotential,  so  that 
Eqc  and  therefore  Ec  were  the  same  for  each  tc.  However,  this  approach  was 
thought  not  to  reflect  the  true  corrosion  state  of  the  alloy  during  the  series 
of  hydrogen  charging/discharging  times  for  a  given  potential.  In  contrast, 
the  method  used  was  considered  to  relate  more  closely  to  the  actual  corrosion 
condition  at  any  time  of  exposure. 

The  drift  in  EQC  (positive  or  negative)  limited  the  data  that  were  valid 
in  terms  of  the  generalized  PDP  conditions  were  somewhat  restricted,  but  the 
flux  control  model  was  found  to  be  applicable  in  all  cases  where  analysis  was 
possible . 

Acetate .  The  anodic  charge  was  again  shown  to  be  almost  entirely 
associated  with  absorbed  hydrogen  (qads  55  10-2  mC  cm-2)  .  Table  15  shows  the 
values  of  J  and  ka  determined  for  charging  times  from  10  to  40  s.  The  value 
of  ka  is  approximately  constant  for  the  three  lower  values  of  T|  but  decreases 
as  Ec  becomes  more  negative.  The  open-circuit  potential  exhibited  a 
significant  negative  shift  at  T|  =  -0.55  V  and  an  even  greater  one  at  -0.6  V 


tv 


m 


m 


(Fig.  10) .  Therefore,  the  decrease  in  ka  is  believed  to  be  an  artifact 
resulting  from  partial  film  reduction/repassivation,  which  inflates  the  value 
of  qa .  The  increased  charge  makes  it  appear  that  less  hydrogen  is  trapped  in 
the  metal  during  egress,  so  there  is  an  apparent  decrease  in  ka .  The  average 
value  of  0.026  s_1  for  ka  taken  over  T}  =  -0.35  to  -0.45  V  is  considered  to  be 
reliable  for  MP35N. 


Table  15 

INTERFACIAL  FLUX  AND  TRAPPING  CONSTANTS  FOR  MP35N 

FOR  ACETATE 


T[  (V) 

_  O  _  I 

J  (nmol  cm  s  ) 

ka  (s'1) 

-0.35 

0.04 

0.029 

o 

o 

1 

0.05 

0.024 

1 

o 

Ln 

0.08 

0.026 

-0.50 

0.08 

0.004 

-0.55 

0.19 

0.0002 

As  for  the  Monel,  the  presence  of  the  film  apparently  results  in  lower 
values  of  J  than  those  for  the  steel.  Values  of  Q'  calculated  using  the 
appropriate  values  of  ka  and  J  for  charging  times  from  0.5  to  40  s  correspond 
well  with  the  experimental  results,  as  shown  for  T)  =  -0.40  V  [Ec  =  -0.425  V 
(SCE) ]  in  Fig.  11.  The  difference  between  the  two  sets  of  data  at  short 
charging  times  is  again  probably  caused  by  changes  in  the  cathodic  current  in 
its  early  stages  and  small  variations  in  EQC  rather  than  a  dependence  of  k  and 
J  on  tc. 


Values  of  the  charge  parameters  for  MP35N  charged  at  T)  =  -0.40  V  are 
shown  in  Table  16.  Both  q.;_  (Fig.  6)  and  q_  (Fig.  9)  follow  the  predicted 
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At  T|  =  -0.65  V,  Eqc  begins  to  have  a  negative  shift  over  the  range  of  tc,  as 
shown  in  Fig.  13,  and  at  -0.7  V,  many  small  bubbles  appeared  on  the  surface. 
These  observations  suggested  that  the  data  for  -0.65  V  were  not  reliable  for 
analysis  in  terms  of  the  hydrogen  ingress  models,  so  the  average  value  of  ka 


was  calculated  excluding  the  result  for  -0.65  V.  The  average  value  of  k&  was 


found  to  be  0.037  s 


Phosphate .  The  variation  in  open-circuit  potential  was  relatively  small 
over  the  range  of  tc  values  (10-40  s)  used  in  the  data  analysis  (Fig.  14) . 
However,  because  qa  generally  varied  in  an  inconsistent  manner  with  tc. 


reliable  values  of  k,  could  not  be  extracted. 

o 


The  values  of  J  and  k_  are  shown  in  Table  19.  The  data  for 

CL 


T|  =  -0.55  V  could  not  be  analyzed  because  qa  decreased  as  t  increased, 
probably  as  a  result  of  the  relatively  large  change  in  Eoc  with  tc  for  this 


overpotential  (Fig.  15).  The  value  of  kg  for  T|  =  -0.65  V  is  undoubtedly 


influenced  by  reduction  of  the  film,  as  reflected  by  the  shift  in  EQC . 
Accordingly,  only  the  results  for  -0.50  and  -0.60  V  were  used  to  determine  the 


average  value  of  ka,  which  was  calculated  to  be  0.058  s 


Table  19 

INTERFACIAL  FLUX  AND  TRAPPING  CONSTANTS  FOR  MP35N 
IN  BICARBONATE 


T\  (V) 

_  9  _  i 

J  (nmol  cm  s  ) 

ka  (s'1) 

1 

O 

cn 

o 

1.08 

0.056 

-0.55 

n 

n 

-0 . 60 

0.49 

0.060 

ID 

O 

1 

0.42 

0.006 

n  =  data  not  analyzable 
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permeation  experiments  J  to  be  ~10  err r  s  .  The  difference  in  diffusivity 

between  this  ESR  steel  and  the  VAR-produced  specimen  used  in  this  study  is 
assumed  to  be  small  enough  to  ignore;  that  is,  it  is  believed  to  be  at  least 
within  the  same  order  of  magnitude.  Data  reported  for  DL  in  pure  iron  show 
considerable  variation,  but  a  value  of  about  5  x  IQ-5  cm2  s-1  at  25°C  appears 

.  o  o  pi 

to  be  reliable.  '  Using  these  values  for  Da  and  D^,  Kr  was  calculated  to 
be  approximately  500,  and  therefore  k  was  4  s-1. 

The  trapping  constant  can  be  regarded  as  the  mean  trapping  frequency  of  a 
single  hydrogen  atom  in  the  metal  lattice  and  can  be  related  to  the  trap 
density  on  the  basis  of  a  simple  model**  involving  spherical  traps  of  radius  d 

O 

and  surface  area  4red  .  This  model  yields  the  following  expression  for  k  in 
terms  of  the  concentration  (N^)  of  traps: 

k  =  4rtd2NiDL/a  (14) 

where  a  is  the  diameter  of  the  metal  atom.  The  diameter  of  the  iron  atom  is 
248  pm  but  a  value  for  d  is  also  required  to  estimate  .  This  in  turn 
presupposes  a  knowledge  of  potential  irreversible  traps.  Such  traps  have  been 
defined^  to  have  a  binding  energy  for  hydrogen  of  >60  kj  mol  ,  and  thermal 
analysis  data15  for  4340  steel  showed  that  MnS  inclusions  in  the  steel  have  a 
binding  energy  of  72  kJ  mol-1. 

The  MnS  inclusions  in  the  test  specimen  were  ellipsoidal,  so  assuming 
that  the  traps  in  the  above  model  are  spherical  is  somewhat  of  an 
approximation.  However,  using  the  value  of  9.2  jjm  for  the  mean  radius  of  the 
MnS  inclusions,  the  density  of  traps  was  calculated  to  be  2  x  10®  m-2 .  In 
view  of  the  assumed  spherical  shape  for  the  traps  and  the  MnS  inclusions,  this 
value  is  in  reasonable  agreement  with  the  actual  concentration  of  inclusions, 

Q  _  O 

2  x  10 3  m  .  The  inclusion  content  should  not  be  altered  by  heat  treatment, 
and  in  fact,  the  value  of  is  the  same  for  both  the  HRC  41  and  HRC  53 
specimens,  further  indicating  that  the  irreversible  traps  determined  from  the 
PDP  measurements  are  stable  entities  such  as  MnS  inclusions. 


J 

t] 
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Other  irreversible  traps  such  as  Fe^C  and  Tie  interfaces  could 
participate,  but  on  the  basis  of  the  value  of  N^,  they  appear  to  have  only  a 
relatively  minor  role.  Although  the  trap  density  was  calculated  using  the 
mean  radius  of  MnS  inclusions,  the  value  of  d  would  have  to  be  unreasonably 
small  (0.01-0.1  pm,  i.e.,  3-4  orders  of  magnitude  smaller  than  the  radius  of 
the  iron  atom)  for  k  to  equal  4  s'1  at  a  value  of  N^  that  approached  the 
typical  levels  of  Fe^C  and  TiC  interfaces  in  steels  (10^4-102^ 


The  average  values  of  ka  for  the  Monel  and  MP35N  are  shown  in  Table  21 
The  values  for  the  Monel  are  independent  of  pH,  so  the  hydrogen  ingress  model 
appears  to  be  valid  irrespective  of  the  electrolyte  provided  that  the  state  of 
the  electrode  surface  (in  this  case,  the  state  of  the  oxide  film)  is  not 
affected  by  either  a  change  of  electrolyte  or  cathodic  charging. 


Table  21 

AVERAGE  VALUES  OF  k&  FOR  MONEL  K500  AND  MP35N 


Electrolyte 

ka 

Monel 

(s'1) 

MP35N 

Acetate 

0 . 020a 

0.027 

h2so4 

b 

0.031 

NaCl/HCl 

0.029 

0.037 

Phosphate 

0.017 

c 

Bicarbonate 

0.019 

0.058 

Averaged  for  AR  and  aged  specimens . 
qa  was  essentially  invariant  with  tc- 
Anomalous  results. 


1 
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The  average  value  of  ka  for  the  Monel  in  the  different  electrolytes  is 
0.021  s-^.  Evaluation  of  the  irreversible  trapping  constant  (k)  from  the 
value  of  ka  requires  diffusivity  data  for  the  pure  Ni-Cu  alloy  to  obtain  DL 
and  also  for  Monel  K500  to  obtain  D,;  the  effect  of  defects  such  as  vacancies 

Cl 

and  edge  dislocations  on  hydrogen  diffusion  in  an  fee  lattice  is  assumed  to  be 
small  enough  to  ignore  (discussed  above) ,  and  therefore  the  minor  elements 
rather  than  defects  in  the  Monel  are  believed  to  determine  the  reversible 
trapping  behavior. 

The  diffusivity  of  hydrogen  determined  from  membrane  permeation 
experiments  is  3  x  10-^®  cm^  s-'*'  in  30  at .  %  Cu-70  at.%  Ni^5  at  25°C  and 
1.48  x  10-^®  cm^  s-^  in  cold  worked  and  aged  Monel  K500  (assumed  to  be  at 
ambient  temperature) .  ^  Although  the  levels  of  Cu  and  Ni  differ  slightly  in 
the  Ni-Cu  alloy  and  the  Monel,  the  error  in  using  the  diffusivity  of  the  Ni-Cu 
alloy  for  DL  is  considered  to  be  negligible.  Accordingly,  using  these  data 
for  Dl  and  Da,  Kr  is  found  to  be  ~2  and  therefore,  the  value  of  k  is 
0.042  s-1. 


Hydrogen  embrittlement  in  Ni-Cu  base  alloys  is  known  to  be  assisted  by 
sulfur  segregated  at  grain  boundaries . ^5  If  the  primary  irreversible  trap  in 
both  Monel  specimens  is  assumed  to  be  sulfur  (considered  as  the  atomic  form  in 
solid  solution  rather  than  as  a  copper-nickel  sulfide) ,  then  the  atomic  radius 
of  sulfur  can  be  used  for  d.  The  value  of  a  is  determined  from  the  weighted 
mean  atomic  diameter  of  Cu  (256  pm)  and  Ni  (250  pm) .  Therefore,  with 
d  =  104  pm,  a  =  252  pm,  DL  =  3  x  10~^  cm^  s_1,  and  k  =  0.042  s-^,  the  density 
of  irreversible  traps  in  the  Monel  (AR  and  HRC  35)  is  found  to  be 

O']  _o  , 

2.6  x  10  m  .  The  actual  concentration  of  sulfur,  calculated  from  the 
density  of  Monel  K500  (8.43  g  cm-5)  and  the  wt%  data  for  sulfur  in  Table  1,  is 

1.6  x  10  S  atoms  m  J  for  the  AR  Monel.  This  value  is  also  assumed  to  be 
valid  for  the  aged  Monel  although  the  distribution  of  sulfur  will  probably  be 
different  due  because  of  segregation  at  grain  boundaries. 

The  difference  of  three  orders  of  magnitude  between  the  value  of  and 
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the  concentration  of  sulfur  atoms  may  indicate  that  the  primary  irreversible 
traps  in  the  Monel  are  not  associated  with  sulfur.  Other  minor  elements,  such 
as  Fe,  C,  and  Mn,  could  well  be  responsible  for  other  types  of  irreversible 
traps.  On  the  other  hand,  the  difference  between  the  experimental  and 
expected  values  of  could  result  from  the  segregation  of  sulfur. 

The  concentration  of  sulfur  segregated  at  grain  boundaries  can  be 
considerably  higher  than  that  in  the  bulk,  so  that  the  traps  at  grain 
boundaries  are  likely  to  consist  of  sulfur  clustered  either  in  the  elemental 
form  or,  more  probably,  as  sulfides.  Hydrogen  as  well  as  sulfur  is  known2®  to 
segregate  to  the  grain  boundaries  in  nickel.  Therefore,  sulfur  clusters  are 
expected  to  be  the  traps  predominantly  encountered  by  the  hydrogen,  so  that 
the  value  of  for  the  Monel  in  fact  probably  reflects  these  clusters.  Thus, 
sulfur,  either  as  the  sulfide  or  in  the  elemental  form,  still  appears  likely 
to  be  the  principal  site  for  irreversible  trapping  in  the  Monel. 

The  surface  film  also  influences  the  diffusion  and  trapping  behavior  of 
the  underlying  alloy.  The  hydrogen  atoms  generated  on  the  surface  of  the 
film-covered  alloy  clearly  must  diffuse  through  the  film  before  entering  the 
alloy.  However,  almost  all  the  diffusion  distance  can  be  shown  to  lie  in  the 
metal  rather  than  in  the  film.  The  mean  distance,  that  the  hydrogen 
diffuses  into  a  metal  (assuming  that  no  interference  with  the  diffusion 
process  occurs)  is  given  by 

X  =  ^(2DLtc)  (15) 

The  diffusion  of  hydrogen  through  the  oxide  as  well  as  through  the  alloy 
during  the  charging  period  is  probably  not  accurately  reflected  by  the 
diffusivity  for  the  alloy  alone.  Nevertheless,  the  diffusivity  of  hydrogen 
through  the  film  is  expected  to  be  low,  and  therefore  could  be  comparable  to 
the  hydrogen  diffusivity  for  the  Monel;  that  is,  they  could  be  within  an  order 
of  magnitude.  Accordingly,  if  DL  =  10-1<^  cm2  s-1,  %  ranges  from  0.45  to 

0.89  |lm  for  charging  times  from  10  to  40  s.  The  thickness  of  the  oxide  film 


anodically  formed  on  nickel  in  aqueous  solutions  is  typically  no  more  than 
1-2  nm,  3  so  if  we  assume  a  maximum  thickness  of  2  nm  for  the  film  on  the 
Monel,  only  0.2%  to  0.4%  of  the  diffusion  distance  involves  the  film  for 
charging  times  from  40  s  down  to  10  s. 

The  distribution  of  sulfur  in  the  oxide  is  likely  to  differ  from  that  in 
the  alloy,  so  that  the  mobile  hydrogen  will  encounter  most  of  the  sulfur  traps 
in  one  part  of  its  overall  diffusion  distance.  Thus,  the  experimentally 
determined  trap  density  will  be  a  weighted  mean  of  the  respective  trap 
densities,  Nm  and  N^,  for  the  metal  and  the  film,  or  more  precisely,  of  the 
respective  number  of  traps  encountered  over  each  part  of  the  diffusion 
distance,  8.  A  weighted  form  of  the  mean  must  be  used  because  of  the 
nonuniform  concentration  profile  of  diffusing  hydrogen.  The  number  of  traps 
encountered  over  the  diffusion  distance  in  a  material,  j,  is  proportional  to 
Nj5j,  so  that  is  assumed  to  be  given  as  a  weighted  mean  by  the  following 
expression : 

Ni  “  twmNm8m  +  wfNf8f]/f(Sm  +  8f )  <wm  +  wf ) ]  (16) 

where  wm  and  wf  are  the  weighting  factors  for  the  metal  and  film, 
respectively.  If  the  weighted  sulfur  concentration  in  the  film  and  therefore 
WfNf  is  low  enough,  the  film  could  contribute  to  the  low  value  of  the 
experimentally-determined  trap  density  (N^)  in  comparison  with  the  expected 
trap  density  (Nm) .  However,  because  the  fraction  of  the  diffusion  distance 
that  involves  the  film  is  very  small,  the  difference  in  sulfur  distributions 
is  probably  a  secondary  effect. 


The  average  values  of  ka  for  MP35N  are  listed  in  Table  21.  Unlike  the 
situation  with  the  Monel,  these  values  change  with  electrolyte.  Some  of  this 
variation  may  result  from  shifts  in  the  open-circuit  potential  over  the  range 
of  charging  times  that  arise  from  changes  in  the  state  of  the  alloy  surface 
due  to  film  formation  (evident  from  the  positive  shift  in  EQC)  and  then  film 
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reduction  (evident  from  the  negative  shift  in  Eoc) •  However,  in  several 
cases,  such  as  1]  =  -0.60  V  in  bicarbonate  and  TJ  =  -0.45  V  in  NaCl/HCl,  the 
open-circuit  potential  was  approximately  constant  with  t_  but  the  values  of  k  = 
were  quite  different  (0.060  s-1  and  0.015  s-1,  respectively).  This  suggests 
that  ka  may  depend  on  the  nature  of  the  electrolyte. 

The  trapping  constant  increases  in  the  following  order:  acetate  <  H2S04  < 
NaCl/HCl  <  bicarbonate.  Except  for  acetate,  this  is  also  the  order  of 
increasing  pH.  However,  as  noted  above,  the  acetate  buffer  is  expected  to 
provide  a  more  constant  supply  of  H+  at  the  electrode  surface  than  do  the 
unbuffered  H2S04  and  acidified  NaCl,  especially  at  the  relatively  low  H+ 
concentrations  used.  Thus,  the  variation  in  ka  with  electrolyte  suggests  that 
kg  depends  strictly  on  the  interfacial  pH,  not  the  bulk  pH.  However,  kg  is 
associated  with  a  process  (trapping)  within  the  metal  and  should  not  vary  with 
interfacial  conditions. 


An  increase  in  pH  would  cause  the  reaction  for  H2  formation  to  change 
from  reduction  of  H+  in  acidic  solutions  to  reduction  of  H20  in  neutral  and 
alkaline  solutions.  However,  although  a  change  in  reaction  may  well  affect 
the  amount  of  adsorbed  hydrogen  and  therefore  the  flux  of  hydrogen  across  the 
surface,  it  would  not  be  expected  to  influence  ka  directly.  The  decrease  in 
H+  concentration  with  an  increase  in  pH  may  have  increased  the  time  needed  for 
the  adsorbed  layer  of  hydrogen  to  reach  its  steady  state  coverage  and 
therefore  the  time  needed  for  the  ingress  flux  to  become  constant.  Such 
behavior  would  violate  the  condition  of  constant  flux  assumed  for  the  hydrogen 
ingress  model  and  may  have  been  responsible  for  the  apparent  change  in  k  ,  but 

a 

if  a  constant  flux  were  not  rapidly  attained  on  MP35N,  the  Monel  should  have 
behaved  similarly  and  shown  an  increase  in  k=;  this  was  not  the  case. 

a 

A  change  in  pH  can  also  affect  the  stability  of  the  surface  film  in  terms 
of  its  thickness  and  composition,  and  this  may  provide  an  explanation  for  the 
change  in  ka .  At  higher  pH  values  up  to  moderate  alkalinity  (pH  9-10),  the 
film  is  likely  to  be  thicker,  so  that  hydrogen  atoms  must  diffuse  through  a 
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greater  distance  before  entering  the  metal  itself.  Thus,  the  thicker  the 
film,  the  more  the  hydrogen  would  be  exposed  to  the  type  of  traps  in  the  film 
over  its  diffusion  path.  This  argument  assumes  that  the  dif fusivities  of 
hydrogen  in  the  film  and  the  alloy  are  similar  on  the  basis  that  the 
diffusivity  for  MP35N  is  low  and  therefore  is  probably  comparable  to  the 
diffusivity  for  the  film,  which  is  also  expected  to  be  low. 


The  mean  trap  density  for  the  metal  and  film  is  given  by  equation  (16) . 

If  it  is  assumed  that  the  dif fusivities  of  hydrogen  in  the  film  and  alloy  are 
similar,  the  measured  irreversible  trapping  constant  can  be  represented  by  the 
following  expression: 


k  =  [wm'km8m  +  wf'kf6f]/[(5m  +  8f)  <wm’  +  wf')] 


where  km  and  kf  are  the  irreversible  trapping  constants  for  the  metal  and 
film,  respectively.  Clearly,  the  measured  irreversible  trapping  constant 


depends  on  the  thickness  of  the  film,  and  therefore  the  film  thickness  could 


account  for  the  change  in  k  with  interfacial  pH.  The  k  for  the  Monel  was 


found  to  be  independent  of  pH,  and  this  is  consistent  with  previous  work 
showing  that  the  thickness  of  the  passivating  film  on  nickel  is  nearly 


independent  of  pH. 4 


Clearly,  this  explanation  is  somewhat  speculative  and  more  information  is 
needed  about  the  nature  of  the  film  on  MP35N  in  the  different  electrolytes  to 
show  that  there  is  in  fact  a  significant  change  in  its  thickness. 

Nevertheless,  if  this  explanation  is  valid,  the  trapping  constant  for  the 
acetate  should  be  least  dependent  on  the  film,  and  therefore,  it  should  be  a 
more  reliable  value  with  which  to  determine  the  density  of  the  principal 
irreversible  traps  in  the  alloy  itself.  Furthermore,  because  the  acetate  has 
a  relatively  low  interfacial  pH,  the  time  required  for  the  flux  to  become 


constant  is  likely  to  be  minimal.  Accordingly,  the  value  of  k,  for  the 

a 


acetate  ha3  been  used  to  calculate  N^,  below. 
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Evaluation  of  k  from  the  value  of  ka  requires  diffusivity  data  for  the 
"pure"  Co-Ni-Cr-Mo  alloy  and  for  MP35N.  However,  these  data  do  not  appear  to 
be  available  for  MP35N,  and  the  most  appropriate  data  for  the  pure  alloy  are 
for  50  at . %  Co-50  at . %  Ni,  for  which  is  5.5  x  10"^  cm2  s-^  at  25°C;23 
because  two  other  primary  elements  are  present  in  MP35N,  there  may  be  a 
significant  difference  in  diffusivity  between  the  Co-Ni  alloy  and  the 
quaternary  base  of  MP35N.  In  the  absence  of  diffusivity  data  for  MP35N,  minor 
elements,  particularly  those  acting  as  reversible  traps,  were  assumed  to  have 
little  effect  on  the  low  diffusivity  of  hydrogen  in  the  "pure"  alloy  lattice 
compared  with  the  effect  in  iron,  where  the  true  lattice  diffusivity  is  five 
orders  of  magnitude  higher.  On  this  basis,  Da  was  assumed  to  be  approximately 
equal  to  DL  and  therefore  k  =  ka. 

The  resistance  of  MP35N  to  hydrogen  embrittlement  has  been  correlated 
with  the  concentration  of  phosphorus  and  sulfur  impurities  segregated  at  its 
major  crystallographic  boundaries.30  If  the  sulfur  and  phosphorus  (in  atomic 
form)  are  considered  the  primary  irreversible  traps,  d  can  be  taken  as  the 
mean  atomic  radius  of  sulfur  (104  pm)  and  phosphorus  (110  pm) .  Likewise,  the 
value  used  for  a  is  the  mean  diameter  of  Co  (250  pm)  and  Ni  (250  pm) . 
Therefore,  with  d  =  107  pm,  a  =  250  pm,  DL  =  5.5  x  10-11  cm2  s_1,  and 
k  =  0.026  s-^,  is  8.2  x  1023-  m-3 . 

The  actual  concentration  of  sulfur  and  phosphorus,  calculated  from  the 
density  of  MP35N  (8.43  g  cm-3)  and  the  wt%  data  for  these  elements  in  Table  1, 
is  8.1  x  1024  S  and  P  atoms  m-3.  The  difference  of  three  orders  of  magnitude 
between  the  value  of  N^  and  the  concentration  of  sulfur  and  phosphorus  atoms, 
as  with  the  Monel,  may  indicate  that  these  elements  are  not  the  primary 
irreversible  traps.  Like  the  Monel,  MP35N  contains  other  minor  elements,  such 
as  Fe,  C,  and  Mn,  as  well  as  S  and  P,  that  could  be  involved  in  irreversible 
traps.  Still,  as  with  the  Monel,  the  difference  is  probably  associated  with 
sulfur,  and  phosphorus,  clusters  at  grain  boundaries  providing  the  main  traps 
encountered  by  the  hydrogen.  In  addition,  the  surface  film  could  contribute 
to  yield  a  mean  trap  density  lower  than  the  expected  trap  density  in  the 
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alloy . 


A  further  factor  in  the  case  of  MP35N  is  the  value  of  DL  used  for  the 
"pure"  quaternary  alloy.  As  noted  above,  for  the  50  at . %  Co-50  at . %  Ni 
alloy  may  well  differ  appreciably  from  the  true  DL  for  a  Co-Ni-Cr-Mo  alloy. 

The  diffusivity  for  46  at .  %  Cu-54  at.%  Ni^  j_s  1.8  x  io-'*-*'1  cm^  s-^,  so  the 
change  from  Cu  to  Co  decreases  DL  by  a  factor  of  2,  and  larger  decreases  are 
likely  with  the  addition  of  Cr  and  Mo.  For  example,  the  addition  of 
20  at.%  Cr  to  Fe  at  27°C  decreases  DL  by  three  orders  of  magnitude although 
the  effect  is  likely  to  be  more  pronounced  for  the  bcc  lattice  of  iron  than 
for  a  close-packed  lattice.  Thus,  it  is  quite  possible  that  the  true  DL  for 
the  quaternary  alloy  is  somewhat  lower  than  the  for  the  Co-Ni  alloy;  if  so, 

a  higher  value  of  would  result.  In  view  of  the  effects  of  the  segregated 

sulfur  and  phosphorus  and  of  a  surface  film,  and  the  lack  of  appropriate 
diffusivity  data  for  the  quaternary  alloy,  it  is  still  likely  that  the  sulfur 
and  phosphorus  impurities  are  the  principal  irreversible  traps  in  MP35N. 

Comparison  of  Ingress  Parameters 

The  flux  of  hydrogen  entering  the  metal  is  somewhat  higher  for  4340  steel 
than  the  Monel  and  MP35N,  largely  because  of  the  presence  of  an  oxide  film  on 
the  two  nickel-based  alloys.  Moreover,  the  value  of  k  is  higher  for  4340 
steel,  even  when  the  rate  constant  approximation  for  MP35N  is  taken  into 
consideration . 

On  the  basis  of  the  values  of  k,  the  principal  irreversible  traps  appear 
to  be  MnS  inclusions  in  the  steel  and,  with  some  uncertainty,  sulfur  and 
phosphorus  in  the  Ni-containing  alloys.  Because  the  traps  in  the  steel  are 
considerably  larger  than  those  in  the  other  two  alloys,  they  should  allow  a 
greater  number  of  hydrogen  atoms  to  bind  per  unit  site,  and  thereby  produce  a 
higher  specific  saturability . 

A  trap  with  both  a  high  binding  energy  and  a  large  specific  saturability 
for  hydrogen  atoms  is  thought  to  be  most  conducive  to  hydrogen 
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embrittlement.  '  In  contrast,  metals  containing  a  high  density  of 
well-distributed  traps  with  a  high  binding  energy  but  a  small  specific 
saturability  for  hydrogen  atoms  should  exhibit  little  susceptibility  to 
interfacial  cracking.  Thus,  the  principal  irreversible  traps  determined  from 
the  trapping  constants  for  4340  steel  and  the  nickel-containing  alloys  are 
associated  with  markedly  different  susceptibilities  to  hydrogen  embrittlement 
The  difference  in  the  traps  and  the  lower  interfacial  flux  for  the  Monel  and 
MP35N  are  probably  key  reasons  these  two  alloys  are  more  resistant  than  4340 
steel  to  hydrogen  embrittlement. 
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CONCLUSIONS 


The  diffusion/trapping  model  for  hydrogen  ingress  was  shown  to  be 
applicable  to  the  steel  and  nickel-based  alloys  provided  that  the  alloy 
surface  remains  essentially  unchanged.  On  the  basis  of  this  model,  hydrogen 
ingress  in  these  alloys  was  found  to  exhibit  the  following  characteristics 
under  the  test  conditions: 

•  The  rate  of  ingress  is  controlled  by  the  flux  across  the  interface  and 
not  by  hydrogen  diffusion  in  the  alloy. 

•  The  rate  constant  for  the  flux  into  4340  steel  is  determined  primarily 
by  electrochemical  characteristics  and  not  by  alloy  microstructure. 

•  Only  about  2-3%  of  hydrogen  atoms  formed  during  charging  enter  4340 
steel  and  Monel  K500  compared  with  20-30%  for  MP35N. 

•  The  proportion  of  hydrogen  being  trapped  increases  with  the  amount  of 
hydrogen  entering  the  alloy. 

•  The  irreversible  trapping  constant  and  the  interfacial  flux  of  hydrogen 
entering  the  alloy  are  higher  for  4340  steel  than  Monel  K500  and  MP35N . 

The  principal  irreversible  traps  identified  on  the  basis  of  the  irreversible 
trapping  constants  appear  to  be  MnS  inclusions  for  4340  steel  and,  with  some 
uncertainty,  sulfur  and  phosphorus  clusters  at  grain  boundaries  in  the 
Ni-containing  alloys.  The  difference  in  both  the  nature  of  the  traps  and  the 
interfacial  flux  are  thought  to  contribute  to  the  difference  in  susceptibility 
of  4340  steel  and  the  two  nickel-containing  alloys  to  hydrogen  embrittlement. 
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Appendix  A 


THE  BASIC  PROGRAM  PDP  173/488 


PAR  1 73/488  is  to  apply  the  potentiostatic  double  pulse  technique 
using  the  PAR  173  Potentiostat  with  276  Interface  module 
Address  of  173/276  =  14 
Terminator = CR  LF  ( &  EOI) 


PRINT  "Set  operating  mode  of  Model  173  to  potentiostat" 


OPEN  "COM1 :9600,N,8,2"  AS  1 
PRINT  #1  ."REMOTE  14- 

DATA -2 
READ  1C 

PRINT  #1  /'OUTPUT  14;  l/E  "1C 

DATA1 
READ  CG 

PRINT  #1  /'OUTPUT  14;  IGAIN  "CG 
DM=10A(-IC)*CG 

DATA  0.01 
READ  OS 

PRINT  #1  /'OUTPUT  14;  MM  2" 

INPUT  "|AE|  (mV)  ";EE 

EE=-EE/1000 

ED=A8S(EE) 

IF  ED<.02  OR  ED=.02  THEN  RR=0:  CR=2.5 
IF  ED>.02  AND  ED<.2  THEN  RR=1 :  CR=25 
IF  ED=.2  THEN  RR=1:  CR=25 
IF  ED>.2  THEN  RR=2:  CR=250 
PRINT  #1  /'OUTPUT  14;  MR  “RR 
EC=EE*10A6/CR 


DIM  ll$(2000),  11(2000),  EC(2000) 

PRINT  "tc  (s)",”Eoc  (V)","  State" 

DATA  0.5,  1, 2.  5,  10,  15,  20,  25,  30,  40 
FOR  1=1  TO  10 
READ  TC(I) 

PRINT  TC(I),; 


/Open  communications  port 
:  'Set  potentiostat  to  REMOTE  with  address  14 


/Set  current  range  code:  10  mA  =  -2 


/Sets  current  measurement  gain 


/Read  offset  from  Ecorr 


/Set  modulation  mode  as  arbitrary  waveform 
/Input  magnitude  of  cathodic  potential  step 


/Set  modulation  range 


/EC  =  voltage  in  counts  based  on  CR  V/count 


/Input  charging  times 


IV 


L«  I 


Mm 


»v 


''A'r 


V./ 


:'Set  time  of  anodic  transient  and  time  base  for  sampling 


IF  TC(I)=.5  THEN  TA=1 :  TB=5000 
IF  TC(I)=1  THEN  TA=2:  TB=5000 
IF  TC(I)=2  THEN  TA=2:  TB=5000 
IF  TC(I)=5  THEN  TA=5:  TB=10000 
IF  TC(I)=10  THEN  TA=10:  TB=20000 
IF  TC(I)=15  THEN  TA=15:  TB=30000 
IF  TC(I)=20  THEN  TA=20:  TB=40000! 

IF  TC(I)>20  THEN  TA=TC(I):  TB=50000! 


N=TC(I)*1000000!/TB  :'Calculate  number  of  cathodic  points  and  anodic  points 

M=TA*1 0000001/TB 

NM=N+M 


PRINT  #1  ."OUTPUT  14;  CELL  O' 
PRINT  #1  ."OUTPUT  14;  READE 
PRINT#1  ."ENTER  14- 
LINE  INPUT  #1,E$ 
EOC=(VAL(E$)/1 000) 

PRINT  EOC,; 

EA=(EOC-OS)*1 000 
PRINT  "Commenced",; 


:'Switch  cell  off 

:'Read  average  open-circuit  potential 


PRINT  #1  ."OUTPUT 
PRINT  #1  ."OUTPUT 

PRINT  #1  ."OUTPUT 
PRINT  #1  ."OUTPUT 
PRINT  #1  ."OUTPUT 

PRINT  #1  ."OUTPUT 


14;  CELL  1" 

/Switch  cell  on 

14;  SETE  "EA 

/Set  potential  EA 

14;  SCV  0" 

/Specify  source  curve  as  0 

14.DCV  2" 

/Specify  destination  curve  as  2 

14;  PCV  2" 

/Specify  processing  curve  as  destination  curve  (2) 

14;  CLEAR" 

/Clear  all  active  points  of  all  curves 

FOR  K=1  TO  N 

EC(K)=EC 

NEXT 

FOR  K=N+1  TO  NM 

EC(K)=0 

NEXT 

PRINT  #1  ."OUTPUT  14.  LC  0,"NM  :'Load  waveform  data  into  Processing  Curve 

FOR  K=1  TO  NM 

PRINT  #1.  "OUTPUT  14;"EC(K) 

NEXT 


PRINT  #1  ."OUTPUT  14;  COPY  2,0' 
PRINT  #1  ."OUTPUT  14;  CLR" 
PRINT  #1  ."OUTPUT  14;  S/P  1" 
PRINT  #1  ."OUTPUT  14;  SIE  1" 
PRINT  #1  ."OUTPUT  14;  TMB  "TB 
PRINT  #1  ."OUTPUT  14;  FP  0” 
PRINT  #1  ."OUTPUT  1 4;  LP  "NM-1 


/Copy  data  to  source  curve 

/Clear  Processing  Curve 

/Specify  1  sample/point 

/Specify  current  to  be  sampled 

/Specify  sampling  time 

/Set  storage  location  of  first  point  as  0 

/Set  storage  location  of  last  point  as  NM-1 


PRINT  #1  ."OUTPUT  14;  NC" 

PRINT  #1  ."OUTPUT  14;  TC" 

PRINT  #1  ."OUTPUT  1 4;  WCD" 

PRINT#1,  "OUTPUT  14;  CELL  0" 

PRINT  #1  ."OUTPUT  14;  DC  0,"NM 

FOR  K=1  TO  NM 

PRINT  #1  ."ENTER  14- 

LINE  INPUT  #1,II$(K) 

ll(K)=VAL(ll$(K))/DM 

NEXT 

NAM$  =  "E"  +  STR$(EE)  +  T  +  STR$(TC{I)) 
OPEN  "O",  #2,  NAM$ 

WRITE  #2,  N;M;NM 
WRITE  #2,  TA;TC(I);TB 
FOR  K=1  TO  NM 
WRITE  #2.  II(K) 

NEXT 
CLOSE  #2 
PRINT  "Executed" 

BEEP 

NEXT 

BEEP:BEEP:BEEP 

END 


initialize  for  curve  acquisition 
:'Commence  curve  aquisition 
:'Wait  until  curve  is  done 
:'Turn  ceil  off 

:'Dump  Processing  Curve 


:*Store  data  on  disk 


WfW'i* 


Appendix  B 


THE  BASIC  PROGRAM  INTERFACE 


REM  INTERFACE  calculates  the  anodic  charge  from  the  current/time  transient 

REM  and  derives  the  best-fit  trapping  constant  and  flux 

REM  for  hydrogen  ingress  under  interface  control. 

REM 

REM 

REM  ERF  approximation  constants 

DATA  0.3275911, 0.254829592,  -0.284496736,  1.421413741,-1.453152027,  1.061405429 
READ  P,  A1,  A2,  A3,  A4,  A5 
DATA  3.1415927,  96487 
READ  PI,  F 

DIM  K(50),  R(50),  SD(50),  QA(40),  11(1600),  IA(1 600) 

PRINT  SPACE$(50)  "SELECT  MODE- 

PRINT  SPACE$(35)  "Calculate  QA  and  QC  for  l/t  transient  . Q" 

PRINT  SPACE$(35)  "Derive  best-fit  trapping  constant  .  K" 

INPUT  "Mode  ";MD$ 

CYCLE:  INPUT  "No  of  to  values  (6)  ";NC 
IF  NC=0  THEN  NC=6 

IF  NC<>6  AND  NC>0  THEN  GOTO  CHARTIME 
DATA  10,15,20,25,30,40 
FOR  J=1  TO  6 
READ  TC(J) 

NEXT 

GOTO  MODE 

CHARTIME:  PRINT:  PRINT  "tc  (s) :" 

FOR  J=1  TO  NC 
INPUT"  ”;TC(J) 

NEXT 

MODE:  IF  MD$="K"  THEN  GOTO  RATE 

REM  Integration  of  Current  T ransient 

REM 

INPUT  "File  Number  ";FF 
INPUT  "Potential  (V)  ";EC 

PRINT  “Code  for  electrode  area  (sq.  cm) :  (1) . 4340  Steel  AR  (1.267)" 

PRINT  SPACE$(41)  "(2) . 4340  Steel  HRc  41/54  (1 .1 1 3)“ 

PRINT  SPACE$(41)  "(3) . MP35N  (1.430)" 

PRINT  SPACE$(41)  "(4) . Monel  K500  AR/HRc  35  (1.267)" 

PRINT  SPACE$(41);:INPUT"  ”;AR 
IF  AR  =  1  OR  AR  =4  THEN  AR  =  1 .267 
IF  AR  =2  THEN  AR  =  1.113 
IF  AR  =3  THEN  AR  =  1.43 


mm 


wwu 


PRINT:  PRINT  "Input  current  offset  for  each  tc:'' 

FOR  J=1  TO  NC 

PRINT  SPACE$(30)  TC(J),;:  INPUT'";  IC(J) 

NEXT 

FOR  L=1  TO  NC 

NAM$="E"+STR$(EC)+"T"+STR$(TC(L))  +  STR$(FF) 

OPEN  T,1,NAM$ 

INPUT  #1,N,  M,  NM 
INPUT  #1,  TA,  TC,  TB 
FOR  K=1  TO  N 
INPUT  #1,II(K) 
ll(K)=-ll(K)/AR 
NEXT 

FOR  K=N+1  TO  NM 
INPUT  #1,II(K) 
ll(K)=(-ll(K)-IC(L))/AR 
NEXT 
CLOSE  #1 

REM  Calculate  cathodic  charge 

QC=0 

FOR  K=1  TO  N 
QC=QC+(II(K)*TB*.000001) 

NEXT 

QC(L)=QC 

REM  Calculate  anodic  charge 

QP=0 

FOR  K=N+1  TO  NM 

QP  =  QP+II(K)*TB*.000001 

NEXT 

QP(L)=QP 

NEXT 

PRINTPRINT  Tc  (S)H  SPC(8)  "Qc  (mC/cmft2)"  SPC(5)  "Qa  (mC/cmA2) 
FOR  L=1  TO  NC 
PRINT  TC(L), 

PRINT  USING  "###.##";QC(L), PRINT  SPC(IO); 

PRINT  USING  "#.###";QP(L) 

NEXT 

BEEP:BEEP:BEEP 

INPUT  "If  charge  analysis  is  required,  enter  A";  CY$ 

IF  CY$  =  "A"  THEN  GOTO  BEGIN 
END 


Analysis  of  Charge  Data 


RATE:  REM 
REM 

IF  K(1)<>0  THEN  GOTO  BEGIN 

PRINT:  PRINT  "Input  anodic  charge  (mC/sq.  cm)  :" 

PRINT  SPACE$(28)  “tc  (s)",  “Qa  " 

FOR  J=1  TO  NC 

PRINT  SPACE$(28)  TC(J),;:  INPUT  "";QA(TC(J)) 

NEXT 

BEGIN:  REM  Determine  trapping  constant  for  minimum  mean  deviation  of  flux 
REM  Initial  trapping  constant  =  0.001  sA-1 

K(1)=.001 
1=1 

SD(0)=100 

TEST1:  IF  CINT(K(I))=1  OR  CINT(K(I))=1 0  THEN  GOTO  PASS  ELSE  GOSUB  CHARG 
IF  SD(I)<  SD(I-I)  OR  SD(I)=SD(I-1)  THEN  K(I+1)=K(I)*10: 1=1+1 :  GOTO  TEST1 
PASS:  KK=K(I-1) 

MT=2 

TEST2:  K(I)=KK*MT 
GOSUB  CHARG 

IF  SD(I)<SD(I-1)  OR  SD(I)=SD(1-1)  THEN  MT=MT+1: 1=1+1:  GOTO  TEST2 

NT=1 

KK=K(I-1) 

TEST3:  K(I)=KK+.1*NT*KK/(MT-1) 

GOSUB  CHARG 

IF  SD(I)<SD(I-1)  OR  SD(I)=SD(I-1)  THEN  NT=NT+1 :  1=1+1 :  GOTO  TEST3 

NT =  1 

KK=K(I-1) 

TEST4:  K(I)=KK-.1*NT*KK/(MT-1) 

IF  K(l)>  0  THEN  GOSUB  CHARG  ELSE  PRINT:PRINT:PRINT’k  <  .0001":  GOTO  FINAL 
IF  K(l)>.0001  THEN  IF  SD(I)<SD(I-1)  OR  SD(I)=SD(I-1)  THEN  NT=NT+1 :  1=1+1 :  GOTO  TEST4 
FINAL:  K(I)=K(I-1) 

GOSUB  CHARG 


PRINT  :PRINT  :  PRINT  SPC(13)  " _ 

PRINT PRINT  SPC(15)  ”tc  (s)  Flux  (nmol/cmA2  s)" 

PRINT  SPC(13) " _ " 

FOR  J=1  TO  NC 
PRINT  SPC(15)  TC(J),; 

PRINT  USING  "  ##.##" ;JJ(J) 

NEXT 

PRINTPRINT  SPC{15)  "Rate  constant  = 

PRINT  USING  "#.####";  K(l) 

PRINT  SPC(  15)  "Mean  flux  =  PRINT  USING  ”##.##”;MV 

PRINT  SPC{15)  "Mean  deviation  =  "SD(I):PRINT:PRINT 
INPUT  "If  refinement  of  rate  constant  is  required,  enter  R  ";RC$ 
IF  RC$  =  "R"  THEN  GOTO  CYCLE 
END 


Calculate  ERF  Functions 


ERF:  REM 
REM 

T=1/(1+P‘Z) 

Q=A2+T*(A3+T*(A4+T*A5)) 

E=1  -EXP(-Z*Z)*T*(A1  +T*Q) 

RETURN 

CHARG:  REM  Calculate  Non-Dimensional  Theoretical  Charge  and  Flux 
REM 

FOR  L=1  TO  NC 
R=K(I)*TC(L) 

Z=SQR(R) 

GOSUB  ERF 

Q1=(1-(.5/R))*E 

Q2=EXP(-R)/(SQR(PI*R)) 

QQ=Z*(1-Q1-Q2) 

JJ(L)=QA(TC(L))/(SQR(TC(L)/K(I))*F*QQ) 

JJ(L)=JJ(L)*1 0000001 
NEXT 

REM  Calculate  Mean  Deviation 
S=0 

FOR  K=1  TO  NC 
S=S+JJ(K) 

NEXT 

MV=S/NC 

SS=0 

FOR  K=1  TO  NC 
SS=SS+ABS(JJ(K)-MV) 

NEXT 

SD(I)=.0001*INT(SS*10000!/NC) 

RETURN 


